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Summary 

A formalism will be presented that allows the transforma- 
tion of two-equation eddy viscosity turbulence models 
into one-equation models. The transformation is based on 
an assumption that is widely accepted over a large range of 
boundary layer flows and that has been shown to actually 
improve predictions when incorporated into two-equation 
models of turbulence. Based on that assumption, a new 
one-equation turbulence model will be derived. The new 
model will be tested in great detail against a previously 
introduced one-equation model and against its parent two- 
equation model. 

Introduction 

Since the emergence of sufficient computational resources 
and adequate computer codes to solve the Reynolds-aver- 
aged Navier-Stokes equations* the turbulence models of 
choice have been either algebraic models like the 
Baldwin-Lomax model (ref. 1), or two-equation eddy vis- 
cosity models. The main shortcoming of the algebraic 
models is the necessity to compute an algebraic length- 
scale which becomes increasingly more difficult as the 
geometry and the flow fields become more complex. Fur- 
thermore, algebraic models fail to account for the impor- 
tant non-equilibrium effects and thereby consistently fail 
to predict the onset and amount of separation in adverse 
pressure gradient flows. The Johnson-King model (ref. 2) 
remedies this deficiency, by introducing a transport equa- 
tion for the turbulent shear stress, but the model still 
depends on the specification of an algebraic length-scale. 

Two-equation models have the advantage that they are 
independent of an algebraic length-scale and can therefore 
be applied to more complex flowfields. The models are 
built on the assumption that at least two scales are needed 
to describe the dynamics of turbulence. A number of dif- 
ferent models has been proposed over the years, but they 
all have in common that a time-scale and a length-scale 
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can be extracted from the two independent variables. 
Among the two-equation models, the k-e model (ref. 3) is 
most widely used in Reynolds-averaged Navier-Stokes 
computations. It solves one transport equation for the tur- 
bulent kinetic energy, k, and a second transport equation 
for the turbulent dissipation rate, e. The model has been 
applied to a large number of flow problems and generally 
produced good results. The main shortcoming of the 
model is its inability to accurately predict adverse pressure 
gradient flows, as was found in the AFOSR Stanford con- 
ference in 1981/82 (ref. 4). Furthermore, the model can be 
numerically stiff when integrated through the sublayer and 
a significant increase of computing time can result from 
the necessity to reduce the time step in order to obtain con- 
verged solutions. 

Recently the monopoly of two-equation models as the 
simplest complete (not depending on the specification of 
an algebraic length-scale (ref. 5)) models has been chal- 
lenged by the re-emergence of one-equation turbulence 
models. While one-equation models have been used ear- 
lier (ref. 6), most of these older models solve an equation 
for the turbulent kinetic energy (or the turbulent shear 
stress) and still depend on the specification of an algebraic 
length-scale in order to calculate the eddy viscosity (see 
however references 7 and 8). The model introduced by 
Baldwin and Barth (ref. 9) solves one transport equation 
for the eddy viscosity and is independent of an algebraic 
length-scale. The Baldwin-Barth model is derived from 
the k-e model and a number of additional simplifying 
assumptions. However, in the course of the transforma- 
tion, several diffusive terms are neglected and one diffu- 
sion term that does not follow from the transformation is 
introduced. The effect of changing these terms can not 
easily be apprehended beforehand, and it turns out that the 
Baldwin-Barth model does perform very differently from 
the underlying k-e model even for simple equilibrium 
flows. The change in the diffusive terms also changes the 
behavior of the model near the edge of shear layers and 
renders the equations ill-conditioned in that region as will 
be shown later. 
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Due to these additional assumptions involved in the deri- 
vation of the Baldwin-Barth model, the connection 
between one- and two-equation turbulence models is not 
entirely clear and the use of one-equation models is often 
associated with a significant loss of generality. 

The aim of the present effort is to establish a firm connec- 
tion between one- and two-equation turbulence models 
that will enable the reader to make a judgement as to when 
the use of a one-equation model is appropriate (or even 
superior to a two-equation model). To achieve this goal, 
the k-e model will be transformed to a one-equation model 
based on a set of clearly defined assumptions. Numerical 
results based on the new model will be compared to those 
of the underlying k-e model and the Baldwin-Barth model. 
The computations will show that the new model performs 
very similar to the k-e model in situations were the under- 
lying assumptions are acceptable. It will also be shown 
that the Baldwin-Barth model gives results very different 
from its parent k-e model. 

It is to be emphasized that the purpose of the present work 
is to establish a clear connection between one- and two- 
equation models of turbulence and not to endorse a new 
model for general use. Due to the close relation to the stan- 
dard k-e model, the new model inherits some of that 
model’s deficiencies. Especially for aerodynamic flows, 
the new model does not perform as well as models specifi- 
cally designed for these applications, like the one-equation 
model of Spalart and Allmaras (refs. 10 and 11) or the 
author’s SST k-0) model (refs. 12-14). It should be pointed 
out that the improved performance of these models is 
achieved by allowing the high Reynolds number form of 
the equations to explicitly depend on the distance from the 
surface. This additional degree of freedom allows for a 
more general calibration which can accommodate both, 
boundary layer as well as free shear layer flows more 
accurately. In contrast, the high Reynolds number form of 
the models tested here is strictly local, a feature that is 
very appealing from a computational standpoint 

This paper has greatly benefited from many important 
comments by Barrett Baldwin, Peter Bradshaw, 
Tom Coakley, and Philippe Spalart I want to especially 
acknowledge Philippe Spalart, who brought to my atten- 
tion the problems of the Baldwin-Barth model near the 
boundary layer edge. 


Turbulence Models 


Transformation of the k-e Model 


This section will present a transformation of the high Rey- 
nolds number version of the k-e model to a one-equation 
model. For simplicity the equations will be written in 
boundary layer coordinates — the general form of the equa- 
tions will be given later. The k-e model reads in boundary 
layer coordinates (x-streamwise coordinate, y normal to 
layer): 
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with the following definition for the eddy viscosity: 
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In order to arrive at a one-equation model, we follow 
Baldwin and Barth and express the time derivative of the 
eddy viscosity by the time derivatives of k and e: 
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Replacing the total derivatives of k and e on the right hand 
side by the right hand sides of equationl 1 gives a single 
transport equation for the eddy viscosity, which, however, 
depends on k and £ as well as on the eddy viscosity: 
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This presents a closure problem with one equation for 
three unknowns. In order to close the equation, two addi- 
tional relations have to be provided. The first one is the 
definition of the eddy viscosity, equation 2, which allows 
one to replace e by the eddy viscosity and the turbulent 
kinetic energy: 
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with c = 0.09 . Note that this relation does not involve 
any assumptions and the resulting equation is still equiva- 
lent to the original k-e model. A second relation is needed 
to eliminate k from the right hand side of equation 4 and 
this relation can obviously not be derived from the k-e 
model (otherwise the k-e equations would be overspeci- 
fied). However, there is a relation readily available, that 
relates the turbulent kinetic energy and the eddy viscosity, 
which has been confirmed for a large number of experi- 
mental boundary layer data and was used by Bradshaw 
et al. (ref. 6) in the derivation of their one-equation model, 
namely (see also ref. 15): 
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where a j is a constant and |-5v| is the turbulent shear 
stress. It has been shown in (refs. 11-14) that Bradshaw’s 
relation (eq. et6), when introduced into the k-(0 model 
leads to a significant improvement of predictions for 
adverse pressure gradient flows. Note that the relation 
between the turbulent shear stress and the turbulent kinetic 
energy that results from standard two-equation models is: 
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using a j = Jc ^ . Equation 7 is not confirmed by experi- 
mental evidence. It is therefore to be expected that the 
introduction of equation 6 will actually lead to improved 
predictions of non-equilibrium flows. 

Since we have a complete set of equations, the one-equa- 
tion model can be derived by straightforward substitution. 
The result is: 


Equation 8 is very complicated and difficult to solve 
numerically. However, the contribution of the terms in the 
last parenthesis of the equation is proportional to the dif- 
ference in the diffusive coefficients of the k- and the e- 
equation. For a number of k-£ models these coefficients 
are equal (ref. 16) and the whole term is exactly zero. The 
importance of this term can be evaluated by changing the 
diffusion constants in the underlying k-e model and inves- 
tigate the influence on the solution. In the logarithmic por- 
tion of a boundary layer, the expression inside the 
parenthesis is exactly equal to zero, so that the investiga- 
tion can be restricted to flows away from the surface. 
Three sets of free shear layer computations have been per- 
formed with different values for the diffusive coefficient 
a £ . Note that the calibration of the standard k-e model for 
a logarithmic layer requires a diffusive coefficient of 
o £ = 1.17 instead of the originally proposed o £ = 1.3 
(with k = 0.41). The results of the computations are given 
in table 1 (o^ = 1 .) 


Table 1. Free shear layer spreading rates for the k-e 
model with different values of o e 


Row 

a e = 1.30 

Q 

II 

o e = 1.00 

Far wake 

0.255 

0.256 

0.257 

Mixing layer 

0.099 

0.100 

0.102 

Plane jet 

0.109 

0.109 

0.110 

Round jet 

0.120 

0.121 

0.124 


It is obvious that changes in the diffusion coefficient a £ 
have only limited influence on the computed results (the 
same is true for a k ). A 30% change in a £ leads to a 
change of only 3% in the spreading rates. The assumption 
that o k = c £ is therefore not very restrictive and the corre- 
sponding terms in equation 8 can be neglected. The second 
assumption in the derivation of the model is therefore: 
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o k = a £ s a (9) 

The resulting high Reynolds number form of the equation 
is: 
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Equation 10 involves the inverse of the v. Karman similar- 
ity length-scale Lvk- 
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The v. Karman length-scale was not very successful when 
used in a mixing length model, mainly because it is singu- 
lar whenever the denominator goes to zero. In the present 
one-equation model, the singularity (of the inverse of Lvk) 
is not a problem, because the destruction term that 
involves Lvk can be limited by any other term that has the 
same dimension, as will be shown later. In the framework 
of eddy viscosity transport models, the v. Karman length 
scale was also utilized in an unpublished model by Bald- 
win (ref. 18), and in a model by Durbin et al. (ref. 19), but 
was never formally derived from the k - 1 model. 


The coefficients of the new model follow directly from the 
k-e model constants: 
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The standard k-e model constants of = 1.44, c e2 = 1.92, 
=0.09==a 1 2 and a K = 1.0 have been used. Note that the 
transformation leads to c 2 = 2/a e = 1.71 with 
a = 1.17 . Since a was chosen to be equal to a k and 
not equal to o £ , the coefficient c 2 had to be slightly re- 
calibrated to match the law of the wall. 


The key to the understanding of the new model lies in the 
comparison of equations 6 and 7. For equilibrium flows 
the two formulations are equivalent and the one-equation 
model will be very close in performance to the k-e model. 
For non-equilibrium adverse pressure gradient flows, 
Bradshaw’s relation, equation 6, is actually better con- 
firmed by experiments than equation 7. For these flows the 
ratio of Production/Dissipation becomes larger than one in 
the outer region of the boundary layer and the k-e model 
will give higher shear stresses than the new model. Since 
the k-e model is well known to overpredict shear stress 
levels for these flows, it is to be expected that the new 
model will lead to improved predictions. For flows with- 
out shear, Bradshaw’s relation, equation 6, has no mean- 
ing, and the new model cannot be expected to give good 
results. An example is isotropically decaying turbulence, 
where one-equation models predict that the eddy viscosity 
stays constant, whereas the k-e model predicts, more real- 
istically, a decay of the turbulent variables. This deficiency 
is associated with the lack of a second scale in the model. 


For shear flows the second scale is provided by the mean 
shear rate. Regions where the mean shear is locally zero 
are bridged by the diffusion and the convection terms. 


Note however, that the overwhelming majority of applica- 
tions of turbulence models is for shear flows, for which 
one-equation models are well suited. 

The main assumption in the derivation of the new model is 
that the turbulent shear stress is proportional to the turbu- 
lent kinetic energy. In standard two equation models this 
assumption is equivalent to: 


Production k = Dissipation). (13) 


used in the derivation of the Baldwin-Barth model. How- 
ever, in that model only the production and dissipation 
terms are transformed based on equation 13. The diffusion 
terms are not transformed exactly. The high Reynolds 
number form of the Baldwin-Barth (BB) model is given 


by: 
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The original Baldwin-Barth model solves an equation for 
the turbulent Reynolds number. In order to allow a one-to- 
one comparison of the constants, it has been reformulated 
here (exactly) as an equation for the eddy viscosity. The 
transformed constants for this model are: 
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They are based on Cjj = 1.2, = 2.0, c^ = 0.09=aj and 

a E = 0.7 for the underlying k-e model. The low Reynolds 
number form of this model can be found in (ref. 9) and is 
not repeated here. 

In order to distinguish the new model, equation 10, from 
the other models in this study, we call it (k-e)iE model 
where the subscript stands for one-equation. The main dif- 
ference between the (k-£)jg and the BB model is the form 
of the destruction term (c 2 term). Although the derivation 
of the BB model starts out from the k-e model, the destruc- 
tion term of that model does not follow from the transfor- 
mation. The link between the BB model and the k-e model 
is broken by the introduction of this term and it is for this 
reason that the two models perform very differently. The 
influence of this term on predicted results will be shown 
later. 
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Low Reynolds Number Terms 


The assumptions leading to the (k-e)ig model are obvi- 
ously not correct in the viscous sublayer, so that the low 
Reynolds number terms of the k-e model cannot be carried 
over to the one-equation model. This is not a great loss, 
because the near wall terms of the k-e model are generally 
complicated and difficult to integrate so that a one-to-one 
transformation is not desirable. 


The purpose of damping functions is to reduce the eddy 
viscosity in the sublayer. In the present model this is 
achieved by reducing the production term as the wall is 
approached and by multiplying the high Reynolds number 
eddy viscosity, by a damping function in order to 
arrive at the corrected eddy viscosity, v f . The damping 
functions are designed in a pragmatic way that ensures 
that the resulting model is numerically stable and does not 
require excessive grid resolution near the surface. The 
present author does not consider it important to match the 
exact asymptotic behavior of the eddy viscosity near the 
surface, since the eddy viscosity is much smaller than the 
molecular viscosity in the immediate vicinity of the sur- 
face and has therefore no influence on the asymptotic form 
of the mean flow variables. Two damping functions are 
introduced, in front of the production term and D 2 into 
the definition of the eddy viscosity: 
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and k = 0.41 . Furthermore the molecular viscosity is 
added into the diffusion term in analogy to the k-e model. 
The coefficient A + is equal to A + =13.5. The complete 
form of the equations is given in the next section. 


General Form of the Equations 


As has been pointe d out by Spalart (ref. 10), alternative 
formulations like JU i -U i . are possible, but for thin 
shear flows there is little difference between them. The 
term involving the inverse of the v. Karman length-scale 
becomes: 


fdndn\ 



v J 


An alternative but numerically more expensive form 
would be: 
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Similarly, all y-derivatives are replaced by their complete 
invariant forms. 


As has been pointed out previously, the inverse of the v. 
Karman length-scale can become singular whenever Q 
goes to zero, leading to an infinite destruction term E k _ . 
In order to prevent this from happening, the destruction 
term is limited by a multiple of the Baldwin-Barth destruc- 
tion term, Egg: 


E u= c 3 E bb tanh 
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with a constant c 3 = 7. Egg is defined as: 
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Equation 23 provides a smooth transition between the two 
formulations whenever goes to infinity. For most of 
the flow E k _ e «c 3 E bb and the original formulation is 
recovered. A less smooth transition could be achieved by 
E le = min(E k _ e ‘, c^E BB ) . The numerical results are 
not sensitive to the constant c 3 . 


In order to arrive at an invariant formulation, all occur- 
rences of the strain rate are replaced by the absolute value 
of the vorticity Q: 
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The final form of the equations is: 
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Numerical Results 

Free Shear Flows 



D = 2jp«(l/^- H )dy (36) 
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Self similar shear layers are very important test cases 
which allow to obtain insight into the performance of tur- 
bulence models, without the need for large computer 
resources. In this section the models will be tested against 
the standard free shear cases, namely a self-similar mixing 
layer, the plane and round jet and the self similar far wake. 
The equations are cast into self-similar form following 
Wilcox (ref. 5), resulting in the following two ordinary 
differential equations for the non-dimensional velocity U 
and the non-dimensional eddy viscosity N: 
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where j=l for the round jet and j=0 for the plane flows. 
The non-dimensional variables are defined as follows: 


Self-similar mixing layer: 


u(x,y) = C/j£/(r|) 

(30) 

v t (x,y) =xU 1 N(t ] ) 

(31) 
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Uj is the velocity of the upper stream (the 
has velocity zero). 

lower stream 

Seif-similar far wake: 



(33) 
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Self-similar jet: 

u(x,y) = (37) 

X v 
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The coefficients in the equations can be obtained from ref- 
erence 5 and equation 3 (Sn=2Si c -S € .) (Note that the coef- 
ficient in reference 5 should be 2U for the round jet): 


Table 2. Coefficients for free shear flows 


Flow 

Su 

S N 

V(T\) 

Far wake 

1/2 

0 

1 

Mixing layer 

0 

-u 

-g£fW)*r 

Plane jet 

1/2U 

-1/2U 


Round jet 

U 

0 



Asymptotic Solution Near Shear Layer Edge 

The analysis of the asymptotic solution near the edge of 
turbulent layer is an important part of turbulence model 
evaluation. It is especially important to determine the sen- 
sitivity of the solution to changes in the freestream values 
specified for the turbulence variables outside the layer. 
Models with solutions that change significantly with the 
freestream conditions are not desirable, because the “cor- 
rect” freestream conditions are not known in most applica- 
tions. The problem has been analyzed in detail in ref- 
erence 17 for two-equation eddy viscosity models. It was 
shown that the k-co model has a severe dependency on the 
values specified for co outside the layer. For example, the 
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maximum eddy viscosity inside a self-similar wake can be 
increased by almost one order of magnitude by reducing 
the freestream value of to. It was also shown that the k-e 
model does not suffer from this ambiguity. There is no the- 
ory available to decide whether a model has a freestream 
sensitivity, but the existence of an asymptotic algebraic 
solution near the boundary layer edge seems to be at least 
a necessary condition for a model to be well conditioned. 

The analysis does not depend on which shear flow is 
selected for the analysis, since the existence of an alge- 
braic solution implies that the terms involving Sy and 
decay faster than the other terms near the edge of the layer 
and can therefore be neglected. The new variable 
f\ = t| - 5 is introduced near the shear layer edge, 8, and 
algebraic solutions of the form: 

dU( y = Ar\ a : N(f\) = Brf (41) 

dr\ 

are introduced into the equations. Straightforward algebra 
shows that the exponents for the (k-e)ig model are: 



Therefore, the velocity and the eddy viscosity approach 
the shear layer edge linearly. It is interesting to note that 
the solution for the k-e model is also linear for 
= c £ = 1 so that the asymptotic behavior of the k-e 
model carries over to the present one-equation model. 

The Baldwin-Barth model does not have a solution of the 
form given by equation 41. However, as pointed out 
before, there is no theory available to show that the exist- 
ence of algebraic solutions is a sufficient, or even a neces- 
sary condition to prevent free stream dependency. 
Numerical test will have to be used to obtain insight into 
the model characteristics. 

Figure 1 shows spreading rates computed with the Bald- 
w in-Barth model as a function of the freestream value, N f 
and the number of grid points, n, across the layer 
0 < t| < 0.4 . The gridpoints are evenly distributed and the 
highest freestream values shown in figure 1 are 2% of the 
maximum value of N inside the layer. The Baldwin-Barth 
model shows a strong sensitivity to the values specified for 
N f , especially as the grid is refined. 

Figure 2 shows the computed velocity profiles on the fin- 
est grid (n=4000) for the highest and the lowest freestream 
values. It is apparent that the model develops extremely 
high gradients in the velocity profile for the lower 
freestream values. The high gradients are the reason why 



Free Stream Value 

Figure 1. Spreading rates of the Baldwin-Barth model 
for far wake, depending on freestream values, Nf, and 
number of gridpoints, n. 

the freestream sensitivity does not show up on the coarse 
grids, where they can simply not be resolved. No grid 
independent solution could be obtained for the low values 
of Nf. Even if the number of points is doubled again to 
n = 8000 , the solution develops even higher slopes and 
lower spreading rates. 

The non-dimensional eddy-viscosities shown in figure 3 
demonstrate even more dramatically the changes due to 
the freestream influence. Like with the k-o model, the 
influence of the freestream values is not confined to the 
immediate vicinity of the boundary layer edge, but can be 
observed throughout the layer. In the present case, the 
maximum of the eddy viscosity changes by about 50% due 
to the changes in the freestream value. Note also that the 
solution for the low freestream values is not grid indepen- 
dent (although 4000 points were used in the computation); 
even lower values of the eddy viscosity can be achieved 
by further increasing the number of points. By contrast, 
the k-co model becomes more diffusive with decreasing 
freestream values and no grid sensitivity is observed. 

Spalart and Ailmaras (ref. 10) have investigated the 
behavior of a turbulent front and found that for the Bald- 
win-Barth model the front propagates into the (physically) 
incorrect direction away from the non-turbulent region. 
They also report that their results are dependent on the 
freestream values specified for the eddy viscosity, consis- 
tent with the present findings. Spalart (ref. 20) also reports 
similar problems with the Baldwin-Barth mode! as shown 
in figurcl for his computation of a self similar mixing 
layer. Like in the present calculations, his results are 
highly sensitive to grid resolution and free stream values. 
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Figure 2 . Velocity profiles for far wake with Baldwin- 
Barth model for highest and lowest ffeestream values 
(number of gridpoints, n=4000.) 



Figure 3. Eddy viscosity profiles for far wake with 
Baldwin-Barth model for highest and lowest 
freestream values (number of gridpoints, n=4000.) 

On the other hand, Baldwin and Barth (ref. 9) have tested 
the model for a flat plate zero pressure gradient boundary 
layer (based on a Navier-Stokes code) and found only a 
moderate dependency on the freestream values. Note how- 
ever, that Navier-Stokes grids are generally fairly coarse 
near the boundary layer edge (assuming that all points are 
plotted in reference 9, there are about forty points across 


the boundary layer but only ten of them in the wake 
region) so that it is very likely that the problem was not 
resolved in that computation. Note also that the present 
author has tested the Baldwin-Barth model in Navier- 
Stokes codes (ref. 21) and did not realize the severity of 
the problem for the same reason. However, Rogers 
(ref. 22) reported that he could not obtain grid independent 
solutions for airfoil flows with the Baldwin-Barth model 
using the 1NS2D Navier-Stokes code. 

Goldberg (refs. 23 and 24) has proposed a number of 
“pointwise” one-equation models, all of which are, in their 
high Reynolds number version, formally identical to the 
Baldwin-Barth model (note that the present problems are 
independent of the low Reynolds number treatment). In 
reference 23 no problems with the model are reported. In 
reference 24 problems in the version given in reference 23 
are conceded and a new version of the model is proposed. 
The author claims that the new model avoids the problems 
near the boundary layer edge due to a regrouping of the 
terms, but fails to explain why the order of the terms in an 
equation would have an impact on the solution. In addi- 
tion, the analysis regarding the behavior of the new model 
near the edge of the layer is incorrect (only an arbitrary 
portion of the coefficient in front of the destruction term is 
considered) and the new model in fact has the same prob- 
lems as the old one. Again the author does not report any 
deficiencies with the new model in the testing portion of 
the paper, although results for self similar shear flows are 
presented. Furthermore, the results presented for the self 
similar flows could not be reproduced with any reasonable 
combination of freestream values and grid distributions by 
the present author. References 23 and 24 are therefore not 
considered in the present discussion. 

What are the implications of the results shown in Figs. 1-3 
for Navier-Stokes applications? There are two different 
strategies. The first one is to specify small values for the 
eddy viscosity in the freestream (inflow). The advantage 
of low values is that they can be specified unambiguously- 
values that are a fraction of the molecular viscosity will 
ensure that they are small compared to those inside the 
layer. However, due to the large gradients developing in 
this case, no grid independent solutions can be obtained, a 
situation not acceptable in a Navier-Stokes code. The sec- 
ond strategy is to specify large freestream values (say x% 
of the maximum values inside the layer). This in turn puts 
a large burden on the user, who has to estimate the levels 
of the eddy viscosity beforehand. Furthermore, the com- 
putations are sensitive to the value of x, so that the results 
depend on an ambiguous quantity. Note also that computa- 
tions over complex geometries involve a number of differ- 
ent turbulent layers, and what is a large freestream value 
for one layer is a small one for an other. Tests have shown 
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that in order to reproduce the experimental spreading rates 
of free shear flows, freestream values of about 20-30% of 
the maximum inside the layer have to be used. Values that 
high would severely impact laminar regions in the flow- 
field and are certainly not acceptable. 

The (k-e)iE model was subjected to the same tests as the 
Baldwin-Barth model. The model did not show any 
freestream dependency, as long as the freestream values 
are small (<1%) compared to the values inside the layer. 
Furthermore the solution for the (k-e) 1E model did follow 
the algebraic solution given by equation 41 near the 
boundary layer edge. Even for very small values, grid 
independent solutions were obtained with only about 
15 points inside the half-layer. 

Table 2 compares the spreading rates of free shear layers 
as computed with the standard k-e and the (k-e) 1E model. 
The experimental values are taken from reference 5. The 
Baldwin-Barth model is not included, because no grid- or 
freestream independent solutions could be obtained. It is 
interesting to note that the (k-e) 1E model gives very simi- 
lar spreading rates to the standard k-e model, despite the 
fact that Bradshaw’s relation, equation 6, is not generally 
true for these flows. The (k-e) 1E model gives somewhat 
lower spreading rates than the standard k-e model for the 
mixing layer and virtually identical results for the far wake 
and the plane jet. For the round jet the new model predicts 
even higher spreading rates than the k-e model so that the 
round jet - plane jet anomaly is somewhat enhanced by the 
new model. 


Equilibrium Boundary Layer Flows 

Wilcox (ref. 5) has popularized the use of defect layer 
computations for testing turbulence models under equilib- 
rium pressure gradient conditions. The equations describ- 
ing these flows are identical to the equations for the free 
shear layers, equations 27 and 28, with the non-dimen- 
sional variables defined as follows: 


u(x,y ) = U e -u x U (•>!); v,(x,y) = N (T|) U f b 



(43) 


with being the friction velocity and 6 the displace- 
ment thickness. The coefficients in equations 27 and 28 
are 


Table 4. Coefficients for defect layer 

s u 

S N 

V(t)) 

p T 

-(l+2P r ) 

-(1 + P T ) -n 


with the non-dimensional pressure gradient parameter 
defined as: 


8 dp 
2 dx 


(44) 


Wall function boundary conditions are specified at the first 
point off the surface (ref. 5): 


Table 3. Spreading rates for free shear flows 


Flow 

k-e 

(k-e) ie 

Experiment. 

Far wake 

0.256 

0.250 

0.365 

Mixing layer 

0.100 

0.084 

0.115 

Plane jet 

0.109 

0.111 

0.100-0.110 

Round jet 

0.120 

0.131 

0.086-0.095 


It is well known that the k-e model gives too low spread- 
ing rates for the far wake and so does the (k-e) 1E model, 
however, the author agrees with reference 8 that the far 
wake analysis should not be given the same weight as the 
other flows, because it is only valid far away from the 
body. The more important near wake is strongly dependent 
on how the wake was generated and is not covered by the 
present analysis. 


C/(n) =-Iln(Ti ) + £/ 0 (45) 

N( T|) = KX\ 

U 0 is obtained from the integral constraint: 

f° U (T|) dr\ = 1 (for details see (ref. 5)). 

Jo 

At the freestream boundary, t/(T|) is set to zero and 
N (T|) to a small freestream value (compared to it’s val- 
ues inside the boundary layer). 

As expected, the Baldwin-Barth model has the same 
freestream dependency as for the free shear layers. Figs. 4 
and 5 show computations with this model for the zero 
pressure gradient boundary layer experiment of Wieghardt 
(ref. 4) on two different grids and with two different 
freestream values. The coarse grid consists of 70 points 
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Figure 4. Velocity profiles for defect layer with 
Baldwin-Barth model and two different freestream 
values, Nf, (fine grid with n=1000). 



Figure 5. Velocity profiles for defect layer with 
Baldwin-Barth model and two different freestream 
values, Nf, (coarse grid with n=70). 

and the fine grid of 1000 points; the high freestream value 
is about 1% of the maximum eddy viscosity value inside 
the boundary layer. Again the solution develops extreme 
gradients near the boundary layer edge for the fine grid 
and the low freestream values. Note that only a moderate 
sensitivity can be observed on the coarse grid. Note also 


that 70 points in the defect part of the boundary layer is a 
very fine grid in Navier-Stokes computations. It is obvious 
that the freestream sensitivity could not be resolved in 
(ref. 9). As in the free shear layer computations, the influ- 
ence of the freestream values and grid densities is not con- 
fined to the vicinity of the boundary layer edge, but affects 
the whole layer. Table 5 shows the computed Cf- values for 
the different cases, reemphasizing this point. 


Tbble 5. Cf results for the Baldwin-Barth model for 
different grids and freestream values, N f 


Gridpoints 

High Nf 

Low Nf 

Experiment 

n=70 

0.0023 

0.0022 

0.0025 

n=1000 

0.0021 

0.0018 

0.0025 


No grid converged solution could be obtained with the 
Baldwin-Barth model for the low freestream value. From 
the present study it is not clear, whether an asymptotic 
solution exists for this model as Nf goes to zero and the 
number of grid points goes to infinity. Whenever the num- 
ber of grid-points was increased, the solution changed 
with an increase in the slope near the edge and a lower 
eddy viscosity inside the layer. 

Figure 6 shows a comparison of results for the zero pres- 
sure gradient case of Wieghardt (ref. 4) for the k-e and the 
(k-e)iE models. Both models give very accurate velocity 



0.0 0.2 0.4 0.6 0.8 1.0 

u/U. 

Figure 6. Velocity profiles for defect layer, Pr=0. 
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profiles and Cf predictions. 

Figure 7 shows results for the adverse pressure gradient 
flow of Clauser (ref. 4) for a non-dimensional pressure 
gradient of = 8.7 . It is well known that the standard 



u/U, 

Figure 7, Velocity profiles for defect layer, pj=8.7. 

k-e model overpredicts the skin friction for adverse pres- 
sure gradient flows, in this case by about 50%. Note that 
the apparent differences in boundary layer thickness 
between the computations and the experiment are a result 
of the definition of T|, involving the friction velocity u^. 
The introduction of Bradshaw’s relation, equation 6, into 
the (k-e)jE model obviously improves the predictions, but 
the skin friction is still too high by about 30%. Due to the 
close relationship of the (k-e) 1E model to the standard k-e 
model, it had to be expected that the deficiency in the 
adverse pressure gradient behavior of the k-e model would 
not entirely be avoided by the new model. 

Navier-Stokes Computations 

All of the following incompressible test cases have been 
computed with the NASA Ames INS2D and INS3D codes 
(ref. 25); the compressible flow was solved with the 
NASA Langley CFL3D code (ref. 26). All flows in this 
study are part of a test base assembled by the author to 
evaluate the performance of turbulence models (refs. 11 - 
14). The flows have been set up in a way to match the 
experimental boundary conditions as closely as possible. 
Furthermore, all computations are performed on grids that 


have been shown to produce grid independent solutions. 

It was initially intended to compare results of all three 
models for the following Navier-Stokes applications. 
However, due to the severe deficiencies discovered in the 
Baldwin-Barth model for the equilibrium flows, the model 
was dropped from the study. It is felt that results obtained 
with that model bear little meaning, since they are invari- 
ably either grid- or ffeestream dependent. 

Flat Plate Zero Pressure Boundary Layer 



Figure 8. Skin-friction coefficient for flat plate 
boundary layer. 

Figure 8 shows a comparison of the computed wall skin 
friction coefficients, Cf versus displacement thickness, 0, 
for a flat plate zero pressure gradient boundary layer. The 
computations are compared with the experimental correla- 
tion of v. Karman and Schoenherr. The low Reynolds 
number version of the k-e model is due to Launder and 
Sharma. The (k-e) 1E model is in good agreement with the 
experiments, whereas the Launder-Sharma model gives 
values about 10% too high. The new model gives accurate 
results, as long as the first grid point satisfies y t + < ~2.5. 

The velocity profiles in inner coordinates are depicted in 
figure 9. Again the (k-e) 1E model is in better agreement 
with the law of the wall than the Launder-Sharma model. 
The good results of the one-equation model are not 
surprising since the low Reynolds number functions have 
been optimized for this flow. Note that the discrepancies 
between the computations of the k-e model and the 
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Figure 9. Velocity profiles for flat plate boundary layer 
in inner coordinates 

experimental correlations are a result of the low Reynolds 
number terms in the Launder-Sharma model. The high 
Reynolds number version of the model with wall functions 
gives the correct skin friction as demonstrated by the 
equilibrium computations shown in figure 6. Results 
similar to the present ones have been reported by Coakley 
(ref. 27) for the Launder-Sharma model. 

Samuel-Joubert Adverse Pressure Gradient Flow 

The Samuel-Joubert flow (ref. 28) is one of the standard 
test cases for non-equilibrium adverse pressure gradient 
flows. In this flow, a flat plate boundary layer develops 
under an increasingly adverse pressure gradient The unit 
Reynolds number is 1.7 x 10 6 . The flow is retarded, but 
not separated. The computations have been performed on 
a 90 x 90 (verified on a 120 x 120) grid. 

Figure 10 shows a comparison of the computed and the 
experimental wall shear stress distributions. The standard 
k-e model gives higher Cf -vaiues than the experiments 
already close to the inflow as a result of its failure to accu- 
rately predict a zero pressure gradient boundary layer 
(fig. 8). The (k-e)i E model is in better agreement with the 
data, but develops somewhat too high wall-shear levels as 
the pressure gradient becomes more severe. This is in 
agreement with the findings for the equilibrium flows. 

Since the wall shear depends on the specifics of the damp- 
ing functions, a better comparison of the high Reynolds 



Figure 10. Skin-friction for Samuel-Joubert flow. 



0 2 

u/U, 

Figure 11. Velocity profiles for Samuel-Joubert flow at 
x=1.16, 1.76, 2.26, 2.87, 3.40 m. 

number behavior of the models can be obtained from the 
velocity profiles. Figure 11 shows that the two models pro- 
duce almost identical results - a testimony to the correct- 
ness of the assumptions leading to the (k-e)i£ model for 
this flow. It was found previously that the Samuel- Joubert 
flow is not a very severe test case and non-equilibrium 
effects are small, so that a good correspondence between 
the two models had to be expected. 
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Driver Separated Adverse Pressure Gradient Flow 


In Driver’s flow (ref. 29), a turbulent boundary layer 
develops in the axial direction of a circular cylinder. A 
strong adverse pressure gradient is imposed on the flow by 
diverging wind tunnel walls plus suction applied at these 
walls. The pressure gradient is strong enough to cause the 
flowfield to separate. The inflow Reynolds number is 
2.8xl0 5 based on the diameter, D, of the cylinder. The 
inflow boundary layer thickness is about 0.2D. The exper- 
iments offer independent wall-skin friction measurements 
and it was found in previous tests (ref. 12) that the data are 
highly self-consistent and well suited to test models under 
strong pressure gradient conditions. The computations are 
performed on a 60 x 3 x 60 (verified on a 90 x 3 x 90) grid. 

Figure 12 shows the wall skin friction coefficient for this 
flow. As in previous comparisons, the standard k-e model 



x/D 

Figure 12. Skin-friction for Driver’s case CSO. 


predicts significantly higher values than the experiment, 
again partly due to the low Reynolds number form of the 
equations. The (k-e)^ model is in significantly better 
agreement with the data, but still somewhat too high, 
especially in the region where the experiment shows 
separation. 

The wall pressure distribution shown in figure 13 is a bet- 
ter indicator of the high Reynolds number performance of 
the models. The one-equation model is a little closer to the 
data, but both models underpredict the viscous-inviscid 
interaction. This can also be inferred from figure 14, 
showing the velocity profiles. Again it can be seen that the 
(k-e) model predicts a stronger retardation due to the 
pressure gradient, but not enough to be in good agreement 



x/D 


Figure 13. Pressure distribution for Driver’s case CSO. 



Figure 14. Velocity profiles for Driver’s case CSO at xj 
D=-0.544, -0.091, 0.363, 1.088, 1.633. 

with the data. The reason for the failure to predict the 
retardation of the velocity profiles lies in an overprediction 
of the turbulent shear stress, as can be seen in figure 15. 
The experimental data are shown in a cartesian system 
aligned with the surface of the cylinder and in a streamline 
coordinate system aligned with the streamline in the mid- 
dle of the boundary layer (u/u max = 0.5 ). It was pointed 
out recently (ref. 30) that a streamline coordinate system is 
the most meaningful system to compare turbulent shear 
stresses; a cartesian system can lead to incorrect conclu- 
sions if the streamlines of the flow are not parallel to the 
surface. In previous computations, it was found very 
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Figure 15. Shear stress profiles for Driver’s case CSO 
at x/D=-0.544, -0.091, 0.363, 1.088, 1.633. 

important to accurately compute the turbulent shear 
stresses ahead of the separation (first two profiles in 
fig. 15) in order to correctly predict this flow. Both models 
give too high values in that region. Note however that the 
(k-£)iE model is in better agreement with the data than the 
two-equation model. As the flow encounters more severe 
non-equilibrium conditions, the differences between Brad- 
shaw’s relation, equation 6, and the relation enforced by 
the two-equation model, equation 7, become more severe 
and the predictions of the two models start to deviate. Fig- 
ure 16 shows the ratio of production versus dissipation as 
predicted by the standard k-E model at the location of the 
maximum turbulent shear stress. This ratio is an indicator 
of the non-equilibrium effects and enters into equation 7. 
Since Bradshaw’s relation is generally more realistic than 
equation 7, the (k-e)^ model gives better results than the 
standard k-E model. 


Backward Facing Step Flow 

The backward facing step is one of the most widely used 
test cases for turbulence model evaluation. While early 
results for this flow indicated that the k-e model underpre- 
dicts the reattachment length by -30%, more recent results 
have shown that the model is off by only about 5%. The 
earlier computations had not enough resolution to accu- 
rately predict the flow. The test case in this study is the 
flow of Driver and Seegmiller (ref. 31). The Reynolds 
number, based on the upstream momentum thickness © is 
Re e = 5, 000 and the ratio of the boundary layer thick- 
ness to step height is about 1.5. The expansion ratio is 
1.125. The computations have been performed on a 
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Figure 16. Ratio of Production/Dissipation for 
Driver’s case CSO as computed from the k-e model. 

120*120 grid with substantial grid refinement near the 
step. The computations are virtually identical to those per- 
formed on a 240 x 240 grid. 

Figure 17 shows the computed and the experimental skin 
friction distributions. The standard k-e model underpre- 



Figure 17. Skin-friction for backward facing step flow. 

diets the reattachment location by about 5% and is gener- 
ally not in good agreement with the data in the separated 
region and near reattachment. Different low Reynolds 
number forms of this model give different skin friction 
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distribution, so that this behavior is certainly a result of the 
low-Reynolds number terms. The (k-e) 1E model is in very 
good agreement with the data. The reattachment location 
is predicted within experimental uncertainty, and there are 
no overshoots near reattachment. Especially impressive is 
the skin friction recovery downstream of reattachment, 
where other models tend to fall more severely below the 
experiments (refs. 11, 12, and 14). 

The velocity profiles depicted in figure 18 show that the 
high Reynolds number differences between the models are 
amazingly small. The velocity profiles are almost identical 



Figure 18. Velocity profiles for backward facing step 
flow at x/h=2, 4, 6.5, 8, 14,32. 

even inside the separation bubble and it appears again that 
the two-equation model does not offer an advantage over 
the one-equation model. Both of the present models fail to 
predict the recovery of the velocity profiles downstream of 
reattachment. This is a general problem with existing 
models and has been obser/ed before (refs. 11, 12, 
and 14). 

Turbulent shear stress profiles are shown in figure 19. 
Again, the experimental profiles are depicted in a cartesian 
and in a shear layer coordinate system. The shear layer 
direction is defined as the direction of the streamline at the 
location where the velocity in the profile is half the differ- 
ence between the minimum and the maximum velocity. 
The results of the two models are very similar. The (k-e) 1E 
model predicts slightly lower shear stress levels due to the 
adverse pressure gradient, again as a results of the intro- 
duction of Bradshaw’s relation, equation 6. It is for this 
reason that the reattachment location is predicted more 
accurately by that model. Both models predict the location 



Figure 19. Shear stress profiles for backward facing 
step flow at x/h=2, 4, 6.5, 8, 14, 32. 

of the maximum shear stress in the recirculation region at 
a different location than given by the experiment — a 
shortcoming also observed with other models. This might 
be partly responsible for the failure of the models to pre- 
dict the correct flow recovery further downstream. 

NACA 4412 Airfoil Flow 

The following flowfield was investigated experimentally 
by Coles and Wadcock (ref. 32). It is a NACA 4412 airfoil 
close to maximum angle of attack at low speed. The chord 
Reynolds number is Re c =1.52 x 10 6 and the angle of 
attack is a = 13.87°. In the experiment the transition was 
fixed at x/c=0.023 and 0. 1 on the upper and lower surface 
respectively. As reported in (ref. 33), a laminar separation 
bubble appears in the calculations on the upper surface, in 
front of the transition point. Laminar separation is not 
observed in the experiments, indicating that transition may 
take place before the trip is reached. The computed results 
are insensitive to the transition location as long as it is 
upstream of the laminar separation point. For the follow- 
ing results the transition location was close to the leading 
edge, as predicted by the models. 

A novelty in the present computations, compared to earlier 
ones (refs, 12 and 33) is that the wind tunnel walls are 
included in the predictions. The computations were per- 
formed on two grids, a 101 x 101 H-grid for the wind tun- 
nel and a 321 x 51 C-grid for the airfoil, linked by the 
chimera capability of the INS2D code. The solid walls of 
the wind tunnel are simulated as slip surfaces. The experi- 
mental data were corrected for the fact that the reference 
pressure p r in the experiments was not equal to due to 
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the pressure tab having been close to the airfoil and there- 
fore not in the undisturbed flow regime. The difference 
between the two pressures is obtained by comparing p r 
and as computed by the code at the location of the 
pressure tab. The difference is largely independent of the 
turbulence model employed. The resulting corrections are 
p~ = p r ( 1-0.06) for the pressure and 
jy„stf r (l + 0.06) for the velocity. 

Figure 20 shows a comparison of computed and 
experimental velocity profiles. The results produced by the 



Figure 20. Velocity profiles for NACA 4412 airfoil, 
0=13.87, Re=1.7xl0 6 . 

one-equation model are better than those of the two- 
equation model for the same reasons as in the previous 
flows. Differences between the models are larger than for 
Drivers’s case, indicating that the flow is even more out of 
equilibrium. 

A significant improvement with the data is achieved com- 
pared to computations without tunnel walls (refs. 12 
and 33), especially with regard to the boundary layer 
thicknesses (which were significantly underpredicted 
before) and the velocities outside the boundary layer. 
However, figure 21 shows that the experimental pressure 
distribution still could not be reproduced accurately. The 
computed suction peak is higher than in the experiment, 
resulting in a more severe pressure gradient along the 
upper surface. The reason for the discrepancies is not 
entirely clear, but is not a result of deficiencies in the tur- 
bulence models, since models which almost duplicate the 
velocity profiles still don’t match the wall pressure. Note 
however, that in the experiment a rectangular test section 
was mounted inside a round wind tunnel, with the walls of 



Figure 21. Pressure distribution for NACA 4412 
airfoil, 0=13.87, Re=1.7xl0. 

the rectangular section extending only to about half a 
chord length upstream of the leading edge of the airfoil. 
The combined blockage effect of the airfoil and of the 
inner test section could have interfered with the flow 
ahead of the airfoil. Since the experimental inflow condi- 
tions into the test section are not known and since uniform 
flow so close to the airfoil is unrealistic, the wind tunnel 
walls simulated in the computations were extended to five 
chords upstream of the leading edge and uniform flow was 
specified. These differences in geometry can very well be 
responsible for the discrepancies in the pressure distribu- 
tions. 

It should however be noted, that the differences in the 
pressure distribution are not responsible for the lack of 
separation in the models, since the computed pressure rise 
is actually higher than the one in the experiments. 


Tiransonic Bump Flow 

The final test case is the transonic bump flow of Bachalo 
and Johnson (ref. 34). In this experiment an axisymmetric 
boundary layer interacts with a shock wave created by a 
circular arc. Only the highest Mach number case 
(Ma=0.925) will be shown. The number of gridpoints was 
150x3x80 which was found to produce grid independent 
results in a previous study (ref. 13). 

Figure 22 shows the wall pressure distribution computed 
by the two models, compared with the experiments. The 
(k-e)iE model gives significantly better results than the 
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Figure 22. Pressure distribution for transonic bump 
flow, M=0.925. 

standard k-e model. Note that this flow is the strongest 
non-equilibrium flow in the study and the differences 
between equations 6 and 7 are therefore largest Again, 
Bradshaw’s relation leads to a significant improvement in 
the comparison, although the shock is still somewhat too 
far downstream with the new model. 

Conclusions 

The connection between one- and two-equation models of 
turbulence has been reexamined. It was found that the 
standard k-e model can be transformed into a one-equation 
model based on only two assumptions. The first one is that 
the diffusion coefficients in the k- and the e-equations are 
the same. By enforcing this condition in the k-e model, it 
was shown that only minor changes resulted from it. The 
second assumption is Bradshaw’s relation that the turbu- 
lent shear stress is proportional to the turbulent kinetic 
energy. This assumption corresponds to 
Production k = Dissipation k in standard two equation 
models and is therefore exact for equilibrium flows. For 
non-equilibrium flows, Bradshaw’s relation is actually bet- 
ter confirmed by experiments than the relation enforced by 
the standard k-e model. The new model was termed 
(k-e) 1E model and tested against the Baldwin-Barth model 
and the standard k-e model. 

Free-shear layer computations have shown that the Bald- 
win-Barth model is ill-conditioned near the boundary 
layer edge. The model does not posses an algebraic solu- 
tion in that region and produces unlimited gradients in the 


velocity as the grid is refined. Furthermore, results are sen- 
sitive to the freestream values specified outside the layer. 
No grid and freestream independent results could be 
obtained with this model and it was for this reason not 
included in the rest of the study. The reason for the failure 
of the model lies in a destruction term that does not follow 
from the transformation of the two-equation model, but 
was introduced on dimensional arguments. 

The findings for the Baldwin-Barth model re-emphasize 
that the behavior of turbulence models near the turbulent- 
non-turbulent interface is one of the most important 
aspects of turbulence modeling. Shortcomings in that area 
are not confined to the immediate vicinity of the interface, 
but change the solution through the whole layer, essen- 
tially nullifying the calibration process. Unfortunately, not 
enough attention is paid to this problem in the derivation 
and calibration of most new models, leaving the door open 
for potentially devastating failures. This is also true for 
Reynolds stress models. 

The new one-equation model does not suffer from these 
deficiencies and shows a very close similarity to it’s parent 
two-equation model near the boundary layer edge. A com- 
parison of the free shear layer results has shown that the 
predictions of the one- and two-equation models are fairly 
close to one another. Both models predict virtually identi- 
cal results for the far wake and the plane jet. Compared to 
the k-e model, the new model predicts about 15% lower 
spreading rates for the plane mixing layer and about 10% 
higher ones for the round jet. 

For zero pressure gradient boundary layers, the new model 
gives almost identical results to the k-e model, except for 
some deficiencies in the low Reynolds number behavior of 
the two-equation model. A number of increasingly stron- 
ger adverse pressure gradient flows has shown that the 
results of the new model improve compared to the k-e 
model predictions as the non-equilibrium effects become 
more important The improved results confirm that the 
main assumption going into the new model is more realis- 
tic than the relation enforced by the k-e model. 

The computations have shown that there is little advantage 
in using a two- over a one-equation model in predicting 
turbulent shear flows. The new one-equation model has 
produced better results than the two-equation model for all 
but the free mixing layer and the round jet. However, due 
to the close relationship to the standard k-£ model, the new 
model still underpredicts the retardation of flows under 
adverse pressure gradient conditions. 


Appendix: Transforming the k-co Model 
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The transformation leading to the (k-e)i£ can be used to 
transform any two-equation model into a one-equation 
model. An example is the k-<o model of Wilcox (ref. 5): 

gM!KKK<*>) 


In the k-(0 model the definition of the eddy viscosity is: 



With the help of equations 3 and 6 the following one-equa- 
tion model can be derived (high Reynolds number form): 
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The constants in this model follow directly from the k-co 
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